Dynamic equivalence between atomic and colloidal liquids 
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We show that the kinetic-theoretical self-difTusion coefficient of an atomic fluid plays the same 
role as the short-time self-diffusion coefficient Ds in a colloidal liquid, in the sense that the dynamic 
properties of the former, at times much longer than the mean free time, and properly scaled with 
Ds, will be indistinguishable from those of a colloidal liquid with the same interaction potential. 
One important consequence of such dynamic equivalence is that the ratio Dl/Ds of the long-time 
to the short-time self-diffusion coefficients must then be the same for both, an atomic and a colloidal 
system characterized by the same inter-particle interactions. This naturally extends to atomic fluids 
a well-known dynamic criterion for freezing of colloidal liquids [Phys. Rev. Lett. 70, 1557 (1993)]. 
We corroborate these predictions by comparing molecular and Brownian dynamics simulations on 
the hard-sphere system and on other soft-sphere model systems, representative of the "hard-sphere" 
dynamic universality class. 



PACS numbers: 82.70.Dd 

One of the fundamental challenges in understanding 
the relationship between dynamic arrest phenomena in 
colloidal systems 1], and the glass transition in simple 
glass- forming atomic liquids Q, is to determine the role 
played by the underlying (Brownian vs. Newtonian) mi- 
croscopic dynamics. It is a widespread notion that col- 
loidal systems constitute a mesoscopic analog of atomic 
systems regarding the relationship between inter-particle 
forces and macroscopic properties 0, |3| • The molecular 
dynamics simulation of an atomic liquid, for example, is 
expected to yield the same equilibrium phase diagram, 
and a similar dynamic arrest scenario, as the Brown- 
ian dynamics simulation of a colloidal liquid, when re- 
ferring to the same model system [1, Q . Important ques- 
tions, however, remain unanswered, even at normal liquid 
states, far from the neighborhood of the conditions for 
dynamic arrest. For example, while it is well-known that 
monodisperse Brownian liquids will freeze when the long- 
time self-diffusion coefficient Dl reaches about O.lxZ^s, 
with Ds being the short-time self-diffusion coefficient 
("Lowen's dynamic freezing criterion" [7]), no analogous 
criterion has been identified for the corresponding atomic 
liquids. 

In the attempt to develop the extension to atomic liq- 
uids, of the self-consistent generalized Langevin equation 
(SCGLE) theory of colloid dynamics [3| , we have discov- 
ered that a well-defined long-time dynamic equivalence 
between atomic and colloidal liquids emerges upon the 
identification of the kinetic-theoretical self-diffusion co- 



efficient of an atomic fluid [9|, as the analog of the 
short-time self-diffusion coefficient Ds of a colloidal liq- 
uid. In this short communication we describe the phys- 
ical foundations of this extended SCGLE theory, which 
also constitute the physical basis of the referred dynamic 
equivalence. One of the most important manifestations of 
the latter is that the ratio D* = Dl /Ds must then be the 
same for an atomic and a colloidal system characterized 
by the same inter-particle interactions, thus naturally ex- 
tending Lowen's freezing criterion to atomic systems. To 
corroborate these predictions we generate and compare 
molecular and Brownian dynamics simulations on hard- 
and soft-sphere model systems representative of the class 
of systems whose long-time dynamics maps onto the dy- 
namics of the hard-sphere fluid, i.e., that pertain to the 
"hard-sphere" dynamic universality class |10|. 

Let us start by considering a model atomic fluid, 
formed by N spherical particles of mass M in a volume 
V, interacting through the pair potential u{r), whose mi- 
croscopic dynamics is described by Newton's equations. 
The fundamental concept upon which we construct our 
theory for these properties is the role of the mean free 
time. To, defined in the kinetic theory of gases as the 
characteristic timescale that provides the crossover from 
the short-time ballistic motion of the atoms to their long- 
time diffusive transport. It is well known that for correla- 
tion times t much shorter than tq, and for distances much 
shorter than the mean free path Iq, all the particles move 
ballistically. For times t much longer than tq, each parti- 
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cle has undergone many collisions, and its motion can be 
represented as a sequence of random (ballistic) flights of 
mean length Iq ^ l/{na'^). As a consequence, for t :s> tq 
each individual particle moves diffusively, with a diffu- 
sion coefficient given by [llj D° — {loY /tq — IqVq ^ 
y^ksT/M/ina^) (where Iq/tq = vq = {kBT/M)^ is the 
thermal velocity). In fact, the rigorous value of D'^ , de- 
termined by the kinetic theory of gases ^ , is 
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This immediately implies that, on the average, the mo- 
tion of each individual particle crosses over from bal- 
listic to diffusive, at the crossover timescale tq, and 
that this crossover behavior must be described by the 
ordinary Langevin equation [ill , [l^ for the instanta- 
neous velocity v(f) of a representative tracer particle, 
Mdw{t)/dt = -C°v(t) + f°(i). In this equation, the 
friction coefficient C° is defined by Einstein's relation, 

= ksT/D^, with 7?° given by the kinetic-theoretical 
result in Eq. ([T]) above. This strongly suggests a deeper 
analogy between the dynamics of Newtonian and Brown- 
ian liquids, and hence, some form of dynamic equivalence 
between atomic and colloidal dynamics, at least regard- 
ing tracer diffusion phenomena. Thus, in both cases the 
relaxation time tq = [M/C"] of the velocity, due to the 
friction force — C°v(<), defines the crossover from ballistic 
(t << To) to diffusive {t » tq) motion. 

The fundamental difference lies, of course, in the phys- 
ical origin of the friction force — C°v(t) and in the defi- 
nition of the friction coefficient C". In a colloidal liquid 
this friction is caused by an "external" material agent, 
namely, the supporting solvent, which also acts as a heat 
reservoir. In an atomic liquid, in contrast, the friction 
force — C°v(t) is not caused by any external agent, but 
by the spontaneous tendency to establish, or restore, 
through molecular collisions, the equipartition of the en- 
ergy available for distribution among the kinetic energy 
degrees of freedom of the system. Thus, in this case 
the friction force — C°v(t) is the dissipative response of 
the system towards the restoration of this partial ther- 
mal equilibrium through collisional heat transport, and 
the spontaneous fluctuations around this equipartition 
is the source of the fluctuating force f°(i). As a conse- 
quence, the corresponding value of is determined by 
the kinetic-theoretical result for in Eq. ([1]), through 
Einstein's relation C° = ksT/D". 

The friction and fluctuating terms of the atomic 
Langevin equation are, of course, a consequence of the 
interatomic collisions of the tracer particle with the rest 
of the particles in the fluid, an effect that could also 
be described by the kinetic terms of the pressure ten- 

sor n (r, t) of the fluid formed by the surrounding par- 

tides. The configurational components of n (r, i), on 
the other hand, correspond to the usual direct force 



-VtE^i w(i"iT) (where "T" stands for "tracer"). This 
force, whose effects must still be taken into account, 
can also be written exactly as / d'^r[Vu{r)]n{r,t), where 
n{r,t) is the local density of particles at position r re- 
ferred to the tracer particle's center [l^. Thus, from 
the superposition of forces, the complete Langevin equa- 
tion for a tracer particle in an atomic liquid actually 
reads Mdv{t)/dt = -Cv{t)+i^{t) + J d^r[Vu{r)]n{r,t) 
By writing the time-evolution equation of n{r,t) 
(linearized around equilibrium), and by substituting the 
solution in this conflgurational contribution of the di- 
rect interactions, one can transform [l^ this direct force 
term, into another friction term plus its corresponding 
fluctuating force, namely, 

= -C°v(i) +f{t)- f dt'ACit - t')v(t') + F(t). 
dt Jq 

(2) 

In this equation the term involving the time-dependent 
friction function A(^(t) describes the mean dissipative 
friction effects due to the collective direct force on the 
tracer particle, whose random component is the Gaus- 
sian stationary stochastic force F{t). The derivation of 
this equation follows step by step the derivation originally 
carried out for colloidal liquids Il3|: it is now explained 
in detail for atomic fluids in Rcf. [15j. One of the main 
products of such derivation is the following approximate 
but general expression for AC* (t) = AC(i) / , in terms of 
the static structure factor S{k) and of the collective and 
self intermediate scattering functions (ISF) F{k,t), and 
Fs{k,t) ^ of the fluid surrounding the tracer particle. 
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which is then also applicable in the present atomic case. 

In order to evaluate AC*(t) we thus need to deter- 
mine F{k,t) and Fs{k,t) at least in the diffusive regime 
t » T°, since the effects described by A(*{t) in Eq. 
^ only manifest themselves at times longer than the 
crossover time For this, one can employ the gen- 
eralized Langevin equation formalism [l3| to derive the 
general memory function expressions for the ISFs of an 
atomic liquid. Such an exercise is described in Ref. jTij . 
As a result one derives general expressions for F{k,t) 
and Fs{k^ t), which in principle describe the dynamics of 
the atomic liquid in the full time-domain, from the bal- 
listic to the diffusive regimes. Requesting consistency 
of such general expressions, with the physical picture 
leading to Eq. 1^ above, these expressions for F{k,t) 
and Fs{k,t) can be written in the long-time (or "over- 
damped" ) limit in terms of the corresponding memory 
functions, C{k,t) and Cs{k,t). The latter may then be 
approximated, as suggested in Ref. [8], by C{k,t) — 
Cs{k,t) = A(fc) A(*{t), so that the resulting long-time 
expressions for F(k, t) and Fs{k, t) can finally be written. 
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in Laplace space, as [16 1 
Fs{k,z)^ 

and 

F{k,z) ^ 
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where the function A(fc) is given by A(fc) — l/[l+(fc/fcc)^], 
with fee = 1.35(27r/cr). 

For a given static structure factor S{k), Eqs. ([3])- 
constitute a closed system of equations, which turns 
out to be identical to the equations that summarize the 
self-consistent generalized Langevin equation (SCGLE) 
theory of colloid dynamics [Sj]. Thus, these equations, 
which constitute an "atomic" extension of the SCGLE 
theory, predict that the long-time dynamic properties of 
an atomic liquid will then coincide with the correspond- 
ing properties of a colloidal system with the same S{k), 
provided that the time is scaled as D'^t and that the 
respective meaning and definition of D'^ is taken into ac- 
count (i.e., that Z?" is given either by the short-time self- 
diffusion coefficient in the case of the colloidal liquid, or 
by Eq. ([1]) in the case of the atomic fluid). One obvious 
example of such long-time properties is the long-time self- 
diffusion coefficient Dl = limt^oo ((Ar(i))^)/6i, whose 
value for an atomic system, scaled as D* = Dl/D^, is 
thus predicted to be indistinguishable from the corre- 
sponding property of the equivalent colloidal system. 

This is an important and distinct prediction of the SC- 
GLE theory, whose accuracy can be readily checked with- 
out even solving the corresponding self-consistent ap- 
proximate equations above, since one can compare molec- 
ular and Brownian dynamics results for the long-time 
self-diffusion coefficient of a given specific system. Thus, 
in Fig. [T] we plot molecular dynamics data for D^^cj)) 
of a hard-sphere fluid both, in the "usual" atomic units 
a{kBT/Mf/^, and scaled as D*{4>) = D l{4>) / D° {4>) , 
with D^{(t)) given by Eq. ([ij. The same figure also 
presents available Brownian dynamics simulation results 
for Dl{(I)) of the hard sphere system without hydrody- 
namic interactions, scaled as D*{(f)) = Dl{(I))/ , with 
D'^ being the ^-independent short-time self-diffusion co- 
efficient of the Brownian particles. Clearly, the "col- 
loidal" and the "atomic" results for D* collapse onto the 
same curve, which we denote by D*jjg{(j)). One imme- 
diate and important consequence of this comparison is 
the extension to atomic liquids of Lowen's dynamic cri- 
terion for freezing: now for both, the atomic and the 
colloidal HS fluids, the condition D'^g{(j)) w 0.1 occurs 
at = = 0.494. 

This long-time colloidal-atomic correspondence is also 
observed in the a-relaxation time Ta{k; cj)), defined by the 
condition Fs{k,Ta) = 1/e, but only in the metastable 




FIG. 1: Long-time self-difFusion coefficient DL{(j}) (main 
figure), and dimensionless a-relaxation time T*(fc) = 
D" {4>)Ta{k- (f>) at ka — 7.1 (inset), of the hard-sphere fluid. 
The results of the molecular dynamics simulations described 
in Ref. [ij for Dl(<^) are expressed in "atomic units" (i.e., 
with d — aiksT /MY^^ , empty diamonds), and scaled as 
D*((^) = DL{(t>)/D°{(j)), with D°{(j)) given by Eq. © (i.e., 
d = D°((^), full diamonds). The corresponding results for 
T*{k) are the full diamonds in the inset. The other full sym- 
bols are the Brownian dynamics simulation results for D* 
from Refs. (triangles) and [l^ (circles), and for r* from 
our present BD simulation results for HS using the method 
described in Ref. ,20i| (squares). The solid lines correspond 
to the solution of the SCGLE theory. 



liquid regime. This is illustrated in the inset of Fig. 
[1] which demonstrates that the molecular dynamics and 
the Brownian dynamics data for T*{k) = k'^D'-Wa{k;<j>) 
evaluated at ka = 7.1 collapse onto a common curve in 
the metastable fluid regime, 0.5 ^ cj), but for (p below 
freezing, both data depart from each other. The rea- 
son for this is that in reality, for 4> ^ 0.5, the dynamics 
of Fs{k,t) at times t « is not yet described by the 
long-time asymptotic expressions in Eqs. (jH) and ([S]); 
instead, Fs{k,t) is closer to its ballistic approximation 
Fsik,t) sa exp[-fc2u2t2/2], so that Ta{k]4>) sa V^/kvo 
(for colloidal fluids, Ta{k\ cj)) w 1/k^D^ at low densities). 

Let us also point out that besides comparing MD 
vs. BD simulation data, one can, of course, also solve 
Eqs. ([3]), (HI), and (O to actually calculate theoretically 
D*fjg{(f) and T*{k). The solid lines in Fig. [T] are the 
result of such numerical calculation, with S{k) provided 
by the Percus-Yevick approximation [ilj with its Verlet- 
Weis correction _22j. In fact, let us mention that the 
data for T*{k) in the inset were employed to calibrate 
the parameter kc in the function A(fc), thus leading to 
the quoted value k^ — 1.35(27r/(T). 

Another manner to express this long-time dynamic 
universality of atomic and colloidal liquids is to plot 
the simulated data of the dimensionless MSD w{t*) = 
< (Ar(t))^ > /61q as a function of the dimensionless 
time t* = t/TQ. From the generalized Langevin equation 
(CLE) in Eq. ^ one can derive the following equation 
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FIG. 2: Main figure: Scaled mean squared displacement 
w{t*) = < (Ar(t))^ > /6lo as a function of the dimen- 
sionless time t* = t/ro, simulated with molecular dynam- 
ics (MD, solid symbols) and Brownian dynamics (BD, empty 
symbols) for fluids with hard-sphere (HS, circles), truncated 
Lennard- Jones (TLJ, squares) and repulsive Yukawa (Y, tri- 
angles) interactions, for states corresponding to (jjHS = 0.1 
and 0.4. Inset: simulation data for the scaled long-time self- 
diffusion coefficient D* of the same systems (and same sym- 
bols) , plotted as a function of the effective HS volume fraction 
<t>HS = 4>Hs{n, T). 



for ■w{t*) 



-t* )w{t* )dt* , (6) 



whose solution satisfies the short- and long-time limits 
w{t* ^ 0) « ^2/2 and w{t* ^ oo) « D*t* , with 
D* = 1/[1 + AC*(z = 0)]. This scaling, however, hides 
the true timescales for atomic and colloidal systems, 
which differ by many orders of magnitude [l5| . For ex- 
ample, the cross-over time tq from ballistic to diffusive 
motion is of the order of one picosecond in an atomic 
liquid and of tens of nanoseconds in a typical colloidal 
system, whereas the crossover time from short to long 
time diffusion, t/ = cP /D'^ (with d — n~^^^), is of the 
order of tens of picoseconds in an atomic liquid and of a 
fraction of a second in a colloidal liquid. Thus, in a col- 
loidal fluid Tj sa 10'' To, and hence, the ballistic short-time 
regime is generally unimportant. Thus, one normally re- 
defines "short times" by first taking the "overdamped" 
limit (i.e., dropping the inertial term in Eqs. ([2|) and ([6])). 
This only changes the short-time limit to w{t*) « t* , but 
leaves unaltered the long-time limit w{t*) « D*t*. Thus, 
except for these short-time differences, the msd of an 
atomic and a colloidal liquid with the same interactions 
and the same S{k) should be indistinguishable at long 
times. This is precisely what is illustrated in Fig. [21 
which plots the simulated data of w(i*) as a function of 
t* for the hard-sphere system at two volume fractions, 
= 0.1 and 0.4. 

Fig. [2] also illustrates the fact that this colloidal- 
atomic dynamic correspondence is not restricted to the 
hard-sphere fluid, but it actually extends over to systems 



with soft repulsive interactions. This is a direct result of 
combining the present colloidal-atomic correspondence 
for the hard sphere system, with another important scal- 
ing rule, which derives from the principle of dynamic 
equivalence between soft-sphere and hard-sphere liquids 
[lot ■ This principle states that the dynamic properties of 
a colloidal liquid at number concentration n and temper- 
ature T, whose particles interact through a repulsive soft- 
sphere pair potential u{r), will be given by those of the 
hard-sphere system with an effective hard-sphere volume 
fraction (j)Hs(n,T) determined by the iso-structurality 
condition, S{kmak;n,T) = Snsikl^ak '^'t'Hs), which re- 
quests that the height of the main peak of the static 
structure factor of the "real" soft-sphere system, and of 
the effective hard-sphere system, coincide. As a con- 
sequence one has, for example, that the msd scales as 
w{t*;n,T) — WHs{t* ]4>Hs{n,T)) and that the curve 
D*jjg {(f) in Fig. [1] becomes a universal curve for all col- 
loidal soft-sphere liquids, provided that its horizontal axis 
refers to the effective volume fraction 0^/5 (n,T). 

The extension of this scaling to atomic systems is 
immediate once the collision diameter a entering in 
the expression for D°{n,T) in Eq. ([T]) is given a 
proper definition for the soft-sphere potential u{r) con- 
sidered. We found that replacing this collision diameter 
by aHs{n,T) = [6(/)//s(n, T)/7rn]^/^ seems to be a sim- 
ple and reasonably universal approximation. To assess 
the accuracy of the resulting general scaling, we gener- 
ated molecular and Brownian dynamics simulation data 
for the truncated Lennard- Jones system, u{r)/kBT = 
T*-i [(cr/r)i2 - 2{(7/rf + l] e[a ~ r), at fixed reduced 
temperature T* = 1, and for the repulsive Yukawa 
potential u{r)/kBT = K exp[~z{r/a — l)]/{r/a), with 
K = 554 and z = 0.149. In both cases, the densities 
were chosen to correspond to effective HS systems at vol- 
ume fractions (pHsin^T) = 0.1 and 0.4. The results for 
w{t*] n, T) in Fig. [2] illustrate that, plotted in this scaled 
manner, the long-time limit of w(t*;n,T) will not dis- 
criminate between atomic or colloidal systems and be- 
tween soft- and hard-sphere interactions. Thus, from any 
of these simulations one should be able to determine D*. 
To illustrate this, in the inset of Fig. [2] we plot the sim- 
ulated data of D*{n,T) for these two soft-sphere model 
systems as a function of 4>hs(P', T), which clearly collapse 
onto the curve D*^g{(f)Hs), represented in the figure by 
the HS MD data. 

In summary, we have shown that at least for model 
liquids whose structure is dominated by (soft- or hard- 
sphere) repulsive interactions, the long-time dynamics of 
atomic liquids is indistinguishable from the dynamics of 
the colloidal systems with the same inter-particle inter- 
actions. As a consequence, just like the equilibrium ther- 
modynamic and structural properties, some dimension- 
less long-time dynamic properties, such as D* and T*{k) 
(the latter only in the supercooled liquid regime), will 
exhibit the same independence from the short-time mi- 
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croscopic dynamics which otherwise distinguishes atomic 
from coUoidal systems. It is reasonable to expect that 
this dynamic universahty will be useful in understand- 
ing, for example, the relationship between dynamic arrest 
phenomena in colloidal systems, and the glass transition 
in simple glass-forming atomic liquids. 
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